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I. INTRODUCTION 



The dynamics of microscopic systems in strong periodic fields forms a problem of funda- 
mental significance, with a vast variety of applications in quantum optics, quantum chem- 
istry, and mesoscopic systems. If the driving field is of a macroscopic nature, for example, 
a continuous-wave laser irradiation, it is appropriate to describe the complete system in a 
mixed quantum-classical way, i.e., to give a full quantum- mechanical account of the central 
system and its energy loss to ambient degrees of freedom (the electromagnetic vacuum or 
weakly coupled internal degrees of freedoms), but to include the field as a classical exter- 
nal driving force. A solution of the dynamics then requires to simultaneously eliminate 
the ambient freedoms and to integrate the equations of motion with an explicit time de- 
pendence. In principle, this can be done exactly using path-integral techniques. However, 
even a partially analytical solution within the path-integral approach is feasible only for the 
very simplest systems in the class addressed, in particular, for periodically driven, damped 
harmonic oscillators or for driven dissipative two- level systems 0. As soon as nonlin- 
ear forces come into play, the path-integral approach requires to resort to extensive and 
sophisticated numerics, such as Monte-Carlo calculations ||, with their own shortcomings. 

In most cases of interest, it is more adequate to make as much use as possible of the 
methods and approximations that have been developed separately for the two problems 
mentioned above, quantum dissipation on the one hand and periodic driving on the other. 
Specifically, it is desirable to combine a Markovian approach to quantum dissipation, leading 
to a master equation for the density operator, with the Floquet formalism that allows to 
treat time-periodic forces of arbitrary strength and frequency. While the Floquet formalism 
amounts essentially to using an optimal representation and is exact Q, the simplification 
brought about by the Markovian description is achieved only on the expense of accuracy. 
Here, a subtle technical difficulty lies in the fact that the truncation of the long-time memory 
introduced by the bath, and the inclusion of the driving, do not commute: As pointed out in 
Ref. ||, the result of the Markov approximation depends on whether it is made with respect 
to the eigenenergy spectrum of the central system without the driving, or with respect to 
the quasienergy spectrum obtained from the Floquet solution of the driven system. In the 
second case it cannot be treated as a system with proper eigenstates and eigenenergies. A 
Markovian approach based on a quasienergy spectrum has been implemented in recent work 
on driven Rydberg atoms || and driven dissipative tunneling [J/J. 

The purpose of the present paper is to investigate these two Markovian approaches to 
damped periodically driven quantum dynamics, with their specific merits and drawbacks, 
for a linear system where an exact path-integral solution is still available: The parametri- 
cally driven, damped harmonic oscillator allows for a very transparent and well-controlled 
introduction of the different approximation schemes at hand. Their quality can here be 
reliably checked since in this system, the quasienergy spectrum is sufficiently different from 
the unperturbed energy spectrum |8j (this feature is in contrast to the additively driven 
harmonic oscillator, where the difference of two quasienergies does not depend on the driv- 
ing parameters ||), and a comparison with the known quantum path- integral solution [|TJ is 
possible. 

Moreover, by switching to a phase-space representation such as the Wigner function, it 
is possible to elucidate the relationship of the quantal results to the corresponding classical 
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Liouville dynamics. Since this relation is particularly close in the case of linear systems, 
this provides an additional consistency check. Therefore, the emphasis of this paper is 
predominantly on the testing and thorough understanding of the available methods. Their 
application to a strongly nonlinear system where analytical path-integral solutions are far 
beyond our present capabilities, will be the subject of forthcoming publications. 

Forming a convenient "laboratory animal" due to its simplicity and linearity, the para- 
metrically driven harmonic oscillator still shows nontrivial behaviour, interesting in its own 
right. We shall give a brief review of the model and its classical dynamics in Section II. 
The central results of the paper, concerning the applicability and quality of the alternative 
Markov approximations, are presented in the course of the quantization of the system with 
dissipation, in Section III. Its last subsection is devoted to a discussion of the asymptotics of 
the quantal solutions, such as the conservative and the high-temperature limits. Section IV 
contains numerical results for a number of characteristic dynamical quantities as obtained 
for the alternative Markovian approaches, and the comparison to the path-integral solution. 
A summary of the various representations and levels of description addressed in the paper, 
with their interrelations, is given in Section V. A number of technical issues are deferred to 
Appendix [A]. Results for an additive time-dependent force in combination with a parametric 
periodic driving are summarized in Appendix [FJ. 



II. THE MODEL AND ITS CLASSICAL DYNAMICS 



For a particle with mass m moving in a harmonic potential with time-dependent fre- 
quency, the Hamiltonian is given by 

*.<*) = + \k(t)x\ (1) 

where k(t) is a symmetric and periodic function with period T. A special case is the Mathieu 
oscillator, where k(t) = m(ujQ + e cos Qt) with Q = 2n/T. Depending on its frequency and 
amplitude, the driving can stabilize or destabilize the undriven oscillation. Fig. |I] shows 
the zones of stable and unstable motion, respectively, for the Mathieu oscillator, in the 
co>Q-£-plane. The equation of motion for a classical particle with velocity-proportional (i.e., 
Ohmic) dissipation in the potential given in ([]]) reads 

1 

x + 7xH k(t)x = 0. (2) 

m 

By substituting x = yexp(— r yt/2), we can formally remove the damping to get an undamped 
equation with a modified potential 

y+(k(t)/m- 1 2 /4)y = 0. (3) 

Already here, on the level of the classical equations of motion, we can apply the Floquet 
theorem for second-order differential equations with time-periodic coefficients. It asserts || 
that Eq. (|3|) has two solutions of the form 

&t)=e i >*<p(t), = &(-*), <p{t + T) = <p(t). (4) 
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The solution &(t) is related to by the time-inversion symmetry inherent in (|3|). Being 
periodic in time, the classical Floquet function ip(t) can be represented as a Fourier series 



ip(t) = c " e 



(5) 



The Floquet index /i depends on the shape of the driving kit) and is defined only modf2. 
There exist driving functions for which [i is complex so that one of the solutions becomes 
unstable (cf. Fig. |1|). In stable regimes /i is real. On the border between a stable and an 
unstable regime, \i becomes a multiple of £1/2 and the solutions £i(t) and £ 2 (£) are not 
linearly independent. For given k(t), the still depend on the damping 7. We denote 
the limit 7 — > of the functions by ^(t). 

The normalization of the c n is chosen such that the Wronskian W, which is a constant 
of motion, is given by 

W = Ut)Ut)-Ut)Ut) = 2i, (6) 

resulting in the sum rule 

00 

E c£(A* + nn) = l. (7) 

71= — OO 

Returning to the original x-coordinate, we find that the fundamental solutions of (0) 
read 

/*(*) = e^ 2 6(t), i = l,2. (8) 

For constant frequency of the oscillator, = const = mug, the Floquet index and the 
periodic function become /1 = (cu 2 — 7 2 /4) 1//2 and <^(i) = (uJq — 7 2 /4) -1 / 2 , respectively, which 
reproduces the results for a damped harmonic oscillator without driving. 

The Green function for Eq. @ is constructed using Eqs. @ and (D), 

G(t,f) =e-^')/ 2 [6(t)6(0 -&(%(*')] /2i (9) 

= e -7(*-f )/2 £ CnCn , gin _ j') + Q( nt _ _ ( 1Q ) 

n,n' 

In terms of this function, the solution of (|2|) with initial conditions x(to) = xq and p(to) = po, 
reads 

u * ^ 9G(t,t ) po 
x(t, t ) = -x — + —G{t, t ). (11) 

Since the potential breaks continuous time-translational invariance, this solution depends 
explicitly on the initial time to- 
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III. THE DISSIPATIVE QUANTUM SYSTEM 



To achieve a microscopic model of dissipation, we couple the system ([!]) bilinearly to a 
bath of non-interacting harmonic oscillators |I| . The total Hamiltonian of system and bath 
is then given by 

H(t)=H s (t) + H SB + H B , (12) 

where 



is the Hamiltonian of N oscillators with masses m u , frequencies u v , momenta p„, and coor- 
dinates x v . The bath interacts with the system via 

N N 2 

H SB = -x £ 9v&v + ^ 2 E (14) 

^ ^ 2m^ 

which couples the system to each bath oscillator v with a strength g,,. The second term in 
Eq. serves to cancel a shift of the potential minimum due to the coupling [|H],|11[]. The 
bath is fully characterized by the spectral density of the coupling energy, 

N 2 

I(u>) = i r y £-£-5(u>-u> v ). (15) 

We choose an initial condition of the Feynman- Vernon type, i.e., at t = to the bath is in 
thermal equilibrium and uncorrelated to the system, i.e. 

P(*o) = Ps(*o) ® PB,eq, (16) 

where pe.eq = exp(— H B / k-oT) / ti-Q exp(— H^/k^T) is the canonical ensemble of the bath and 
k&T Boltzmann's constant times temperature. 



A. Interaction picture and perturbation theory 

Due to the bilinearity of the system-bath coupling, one can always eliminate the bath 
variables to get an exact, closed integro-differential equation for the reduced density matrix 
ps = trap, which describes the dynamics of the central system, subject to dissipation fl2|-|l4 



In most cases, however, this equation cannot be solved exactly. In the limit of weak coupling, 

7 < k B T/h, (17) 
7 < Aa/3, (18) 

it is possible to truncate the time-dependent perturbation expansion in the system-bath 
interaction after the second-order term. The quantity 7 denotes the effective damping of 
the dissipative system, and A a p are the transition frequencies of the central system (see, 
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e.g. Eq. (|56|), below). The autocorrelations of the bath decay on a time scale h/k^T, and 
thus in this limit instantaneously on the time scale I/7 of the system correlations. 

With the initial preparation (Jl6|), the equation of motion for the reduced density matrix 
in this approximation is given by (T^-|l4| 



H s (t),ps(t) 



h 



tr B 



HsB,ps(t) 



-WJo dTtTB 



SB, 



The tilde denotes the interaction picture defined by 

o(tj) = ui(t,if)6u (t,tr), 

U (t, t') = T exp (-^ £ dt"(H s (f) + H B ] 
where T is Wick's time-ordering operator. 



For H$ and Hsb as in Eqs. (|13|), ([14]), we find the master equation 

fa(t)=-~[H s {t),Ps(t\ 

1 N f°° r 
-tfHd / dr[S u (T)[x, [x(t-r,t),p s (t)]] 

+iA„{r)[x, [x(t-T,t), Ps {t)} + ]}, 

with [A, B]+ = AB + BA and 



A v (t) 



coth ( — 1. I coscu^t, 



2m v uj v 



2k B T 



2m u u! u 



sin a;^, 



(19) 

(20) 
(21) 



(22) 

(23) 
(24) 



the symmetrically ordered and antisymmetrically ordered, respectively, correlation functions 
of the bath oscillator v. 



B. Markov approximation with respect to the unperturbed spectrum 

So far, we have followed the standard approach to dissipative quantum dynamics in the 
weak coupling limit []12|-[14|]. In the following subsections, we shall contrast a simpler Markov 
approximation based on the unperturbed spectrum, with a more sophisticated approach that 
accounts for the modification of the spectrum due to the driving. 



1. Master equation 



In the following, we restrict ourselves to an ohmic bath, 
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I(uj) = mjuj, 



(25) 



fixing the relation between the macroscopic damping constant 7 and the microscopic coupling 
constants g n introduced in Eq. flU]). By imposing a Drude cutoff I(uj) — > — iuj/ud) 

with u D ^> u , Q, divergent integrals are avoided. 

In a crudest approximation, the time dependence of the system Hamiltonian is neglected 
in the derivation of the master equation, i.e., the incoherent terms in the master equation 
are calculated replacing Hs(t) by H$ = 0-/T) J^dtHsit), i.e. the Hamiltonian with zero 
driving amplitude. The position operator in the interaction picture is then given by 

~ P 

x(t, if) = x cosu (t — if) H sinu; (t — t'). (26) 

Since the information on the phase of the driving is lost, it depends only on the difference 
t — if of its arguments. 

Inserting this operator and the correlation functions (|23"D , (pjj) into Eq. (0), leads to the 
master equation 

Ps = ~[H 8 (t),ps] - ^-7 [x, \p,Ps]+] 

--^D pp [x, [x,ps]] + T2D xp [x, [p,p s ]]- (27) 

The right-hand side of this equation depends on time only through its first, the Hamiltonian, 
term and therefore retains the periodicity of the system Hamiltonian Hs(t). 

This form of the master equation does not produce a positive semidefinite diffusion ma- 
trix. It consequently does not exhibit Lindblad form |TB,IB-|r7|]. The positivity of ps is 



thus not guaranteed for all elements of the function space of density operators. The Marko- 



vian approximation implies that quantum effects on a length scale I < X^b = h/ ' \jAmk^T 
(non-Markov effects) cannot be described selfconsistently [fT7-l9||. Note also that within a 



Markov approximation, the master equation is periodic with the driving period T = 2it/£l 
(Floquet form). This is in contrast to the non-Markovian exact master equation []]]. In this 



latter case, the effective master equation has the structure of (p7|) with time-dependent co- 
efficients D xp and D pp that depend also in a non-periodic way on the time elapsed since the 
preparation at t . In Wigner representation, this corresponds to a time- dependent diffusion 
coefficient (see below). 

The coefficients 7 and D pp can be evaluated straightforwardly |2(J to give 



7 = 7, (28) 

-mhun coth - ° . 
2 u 2k B T 



D pp = 7^ m ^o coth ■ Kr %. (29) 



The evaluation of the cross-diffusion D xp is more complex. Because we did not find it 
in the literature, we give the outline of its derivation. The logarithmic divergence of D xp is 
regularized by the Drude cutoff to obtain 




U» = - 5= P/ d.coth| — ]^_ ! _^_, (30) 
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where P denotes Cau chy's principal part. The integral in Eq. (j30D is solved by contour 
integration in the upper half plane. Expressing the resulting sums by the psi function 
ip(x) = dlnT(a;)/da; |2T| and neglecting terms of the order uj /ujd, we obtain 



D 



xp 



h 

71 



i> 1 + 



27ik B T 



+ C 



(31) 



where C is the Euler constant. 

Interestingly enough, m^D xp coincides with the Drude regularized divergent part of the 



stationary momentum variance of a dissipative harmonic oscillator |22 



It must be stressed that the dissipative terms in the master equation ( p7|) are independent 
of the driving. This manifestly reflects that the time dependence of H§(t) has not been taken 
into account in the incoherent terms of the master equation. 



2. Wigner representation and Fokker-Planck equation 

In order to achieve a description close to the classical phase-space dynamics, we discuss 
the time evolution of the density operator in Wigner representation. It is defined by |23 



1 r 00 

W(x,p,t) = — dx'e 2ipx/h (x-x'\p s (t)\x + x'). (32) 
wn J-oo 

The moments of the Wigner function are the symmetrically-ordered expectation values of 
the density operator. 

Applying this transformation to the master equation (p2[), we obtain a c-number equation 
of motion, 

d t W(x,p,t) = L(t)W(x,p,t), (33) 

with the differential operator 

L(t) = — -pd x + ~fd p p + k(t)xd p + ^D pp d 2 p + jD xp d x d p . (34) 

777- 

Equation ( [34]) has the structure of an effective Fokker-Planck operator. However, for D xp ^ 
0, the diffusion matrix is not positive semidefinite; correspondingly ( |33l) has no equivalent 
Langevin representation. 

As is the case for the master equation from which it has been derived, the coefficients 
of the Fokker-Planck operator retain the it periodicity of the driving, so that Eq. ( p3[ ) has 
solutions of Floquet form. This fact will be exploited in the following subsection to construct 
the solutions. 



3. Wigner- Floquet solutions 

The Fokker-Planck equation for the density operator in Wigner representation, Eq. fl33|) 
with Eq. (j34j), offers the opportunity to make full use of the well-known and intuitive results 
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for the corresponding classical stochastic system. In particular, a solution of the Fokker- 



Planck equation can be obtained directly by solving the equivalent Langevin equation fl24| , |25 



or by using the formula for the conditional probability of a Gauss process p5 |. In the present 
case, however, the fact that the diffusion matrix of ([34]) is not positive semidefmite requires 
to take a different route. 

Since Eq. ([33]) with Eq. (0) represents a differential equation with time-periodic coeffi- 
cients, it complies with the conditions of the Floquet theorem. Consequently, there exists a 
complete set of solutions of the form 



W a (x,p,t) = e Ata *M a (x,p,t), u a (x,p,t) = u a (x,p,t + T), 



(35) 



henceforth referred to as Wigner- Floquet functions. 

We construct a solution for (|33"D of this form with fi Q0 = by the method of characteristics 
26| , cf. the Appendix In the limit to ~* — oo, the terms in the first line of ( |A18|) , which 



contain the initial condition, vanish and we obtain the asymptotic solution 



W 00 (x,p,t) 
with the variances 



1 

2^ 



0~xx 

(t) a xp (t) 



-1/2 



exp 



x 
P 



(t) a xp {t) 
O'xpit) o-pp{t) 




o~ xx (t) 
o~ xp (t) 

a pp (t) 



2 1 a 



I'p 



m 

2 1 a 



dt' [G(t,t')} 2 
d 



m 



J^dt'G(t,t')-G(t,t') 



rwyDxp + 2 , ~fD 



pi* 



dt' 



(36) 

(37) 
(38) 

(39) 



Note that in 



the difference in using D pp and D = D pp + 



'xp 



sec 



in 



Appendix |A|) is meaningless, since it is a correction of order 7. By inserting for G(t, t') the 
Fourier representation (0), one finds that the variances are asymptotically time-periodic. 

Starting from Wqo, we construct further Wigner- Floquet functions: By solving the char- 
acteristic equations (see Appendix |A|), we find the two time-dependent differential operators 



Qi+{t) =h{t)d x + mf x {t)d p , 
Q.2 + {t)=f2{t)d x + mf 2 {t)d p . 

They have the properties 

[Lit) - d t , Q 1+ (t)) = \L{t) - d t , Q 2+ (t)) = 

and 

Q 1+ (t + T)=e^ 2+i ^ T Q 1+ (t), 



Q 2+ {t + T)=e 



(-7/2-i/u 



)T Q2 + {t). 



(40) 
(41) 



(42) 



(43) 
(44) 



Taking the commutation relation ( f4"2"D into account, the functions 
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W nn ,(x,p,t) = Q 7 { + (t)Q%' + (t)Woo(x,p,t), n,ri = 0,1,2,... (45) 



also solve Eq. (p3|). 

Due to Eqs. (|43l) , (f44"l), they are of Floquet structure with the Floquet spectrum 

AW = n(-7/2 + i/i) + n'(- 7 /2 - i/i). (46) 

This spectrum is independent of the diffusion constants, as expected for an operator of type 
f2?|| , and therefore is the same as in the case of classical parametrically driven Brownian 



oscillator 

The expression for the eigenf unctions in the high-temperature limit of the (undriven) 
classical Brownian harmonic oscillator in Refs. p7|p9| is also of the structure (f|5|). We can 



recover this solution by inserting the classical diffusion constant mk^T and the undriven 
limit e — > for the classical solution, given in Sect. ||. 

C. Markov approximation with respect to the quasienergy spectrum 

The master equation ( p7|) can be improved by including the time-dependent term in the 
system Hamiltonian (fj) before a Markov approximation is introduced, to account for the 
change in the quasienergy spectrum due to the driving. 

1. Floquet theory and quasienergy spectrum 

For a Schrodinger equation with time-periodic system Hamiltonian such as (|l|), the Flo- 
quet theorem M asserts that there exists a complete set of solutions of the form 

\Mt)) = e-^lMt)), \Mt + T)> = \Mt)) (47) 

The quasienergy \i a plays the role of a phase and therefore is only defined modf2, cf. 
Ref. ||. We shall use the basis {\4> a (t))} as an optimal representation to decompose states 
and operators. 

For the parametrically driven harmonic oscillator ([I]), the Floquet solutions for the 
Schrodinger equation are derived in the literature in various ways |30|-p3|. We skip the 
derivation and merely present the result, 

Mx > t] = [0m) (f§) ' Ha ( x ^ m,nmm ) exp (#(*)*7 2 tf w) » ( 4§ ) 

for the Floquet solutions in the stable regime, where H a is the a-th Hermite polynomial, 
a = 0, 1,2, . . .. The Floquet index for this solution is fi a = fi(a + 1/2). This gives the 
quasienergy spectrum 

H*,k = (a + l/2)/i° + fcfi, k = 0, ±1, ±2, (49) 

Note that ( ^5| ) are solutions only in the stable regime. Consequently ji is real, cf. Sect. [TT]. 
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In analogy to the annihilation and creation operators for the undriven harmonic oscilla- 
tor, one can define operators T and F> which act as shift operators for the Floquet states, 
i.e. 



t(t)\Mt)) = V^\^a-l(t)), 
ft(t)^ a (t)> = V^ + T|Va+l(*)). 



(50) 
(51) 



For a parametrically driven harmonic oscillator, T(t) can be expressed in terms of position 
and momentum operator as [plj|32 



fw=il*VxS(«)-*i 



mh 



(52) 



The relations (pOf) and ([H]) can be proven by inserting the Floquet solutions (j48|) and using 
the recursion relations for Hermite polynomials f2~I |. 



The matrix element X a p(t) of the position operator x with the states \ip a (t)), which we 
shall need later, reads 



with the transition frequencies 



k 

' fT dte~ iknt (Mt)\x\Mt)}, 



T Jo 



A 



a/3k 



(53) 
(54) 

(55) 



(56) 



For Eqs. (|54| ) and (p5[), the periodicity of the Floquet states \4> a (t)) has been used. The 
Fourier components X a pk are preferably evaluated in the spatial representation, 



X a p(t) 



dx ijj a (x, t) xipp(x, t) 



(y/P&WaJI-l + V^l(*)<W 



(57) 
(58) 



by inserting the Fourier expansion ([5]) for £,f(t), to give 



(59) 



2. Improved master equation 



We start anew from the full master equation in the weak-coupling limit, 
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P=~[H a (t),P] 



+— / duj I(u)n t Joj) / dre iWT [x(t-r,t)p,x] + h.c. . (60) 

Tift J-oo Jo 

Here, h.c. denotes the hermitian conjugate of the dissipative part and 

n th (uj) = (e^ T - l)" 1 = -n th (-cu) - 1 (61) 



gives the thermal occupation of the bath oscillator with frequency u. To achieve a more 
compact notation, we have required that = which for an Ohmic bath, cf. 

Eq. (|25|), is just the analytic continuation. 

The fact that the Floquet states \ip a (f)) of the undamped central system, Eq. ([471). 
solve the Schrodinger equation, allows for a substantial formal simplification of the master 
equation: With the density operator being represented in this basis, 

p a p(t) = (Mt)W)\Mt)), (62) 
the master equation takes the form 
1 f°° 

Pa/3 = — r / du I{u)n t h(v) 
TinJ-oo 

x f°°dT^ T £ {X aa r{t - T) Pal0 ,Xy{t) - X*, a (t)X a//3 ,(t - T) Pf3/f3 } + h.c. (63) 

J ° a' (3' 

Inserting fl55|) and (59) and using the identity / °° dr e lulT = 7t5(u>) + P(i/u), we arrive at 
the explicit equation of motion 

^ = JEE{- /(A Q ,^ fc 0^h(Aa'/3'fc0e i(AQ '^ fc '" AaQ ' fc) *^ afe X Q / / 3/^p /3 /^ 

11 a'(3' kk' 

+/(A^ fc )n th (A aQ , fc )e i(A -^- A ^''=' ) 'X Q , aA; p a ^/X; /3 , fc ,} + h.c. (64) 



The quasienergies of the undamped central system appear in Eq. fl64|) by way of the Aq^. 
Since these frequencies contain only differences of quasienergies, they have a direct physical 
significance as transition frequencies and so may be used as arguments of I{u) and nth(<^)- 
This is not the case for the quasienergies themselves, due to their Brillouin-zone-like ambi- 
guity, cf. Eq. P9|). Shifts of the A a/3k brought about by the principal parts of the integrals 
have been neglected. 



3. Rotating-wave approximation and solution in the Floquet representation 

In a rotating-wave approximation (RWA), it is assumed that phase factors exp[i(A Q( g& — 
A a 'f3>k')t], with (a,/3,k) 7^ (a',/3',k') in Eq. (|4j) oscillate faster than all other time depen- 
dences and hence can be neglected. This argument applies, however, only to quasienergy 
spectra without systematic degeneracies or quasidegeneracies. Indeed, the harmonic poten- 
tial we are presently dealing with has the peculiarity of equidistant (quasi-) energy levels, cf. 
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Eg. (f49|) , so that additional terms have to be kept. Here, the condition (a— (3, k) = (a'—/?', k') 
is sufficient to ensure A a p k = A a /^ fc /. Therefore these terms have to be kept in RWA. 

Making the RWA, substituting Eq. (|59|) in Eq. (j64|) , and assuming an Ohmic bath as 
above, we obtain the time-independent master equation 

Pa/3 = 1{( N + !) ^vWlX/^ + lW./m - (a + P)Pafl) 

+N (2fitfp a -. liP -. 1 -{a + (3 + 2)p a ^ | . (65) 

The effective thermal-bath occupation number 

N = E( C °) 2 ' (p° + kQ)n th (p° + kQ) (66) 

k 

reduces to N = ri t h(^o) in the undriven limit. 

Formally, this master equation coincides with that for the undriven dissipative harmonic 
oscillator in rotating-wave approximation [Pfl. It has the stationary solution 

« = IvTT GvTLW (67) 

The density operator of the asymptotic solution is diagonal in this representation and 
reads 

oo 

pUt) = Y,p«a\Mt))(Mt)\- (68) 

The basis {\?pa(t))} corresponds to the "generalized Floquet states" introduced in Ref. ||, 
i.e., they are centered on the classical asymptotic solution and diagonalize the asymptotic 
density operator. 

To get the variances of fl58p, we switch to the Wigner representation, 



WUx, P ,t) = J2pZW a (x,p,t), (69) 



a=0 

where 



W a (x,p,t) = y -^e~ z L a (2z 2 ), (70) 

7T 



z n 



is the Wigner function corresponding to |^ a (t)) [j33|| , with the Laguerre polynomial L c 
Using the sum rule [[Ti] 



£ K a L a (x) = (1 - k)- 1 exp ( J^Lj , (72) 



a=0 

we obtain the asymptotic solution in Wigner representation as 
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(73) 



It is a Gaussian with the variances 



a xp (t) = h(N + 1/2) (Q(t)${t) + (t)g(f)) /2, 
( x pp (t) = /im(iV + l/2)^(t)4 W- 



(74) 

(75) 
(76) 



To enable a comparison between the different equations of motions for the dissipative 
quantum system, we give for the master equation in RWA fl65|) also the corresponding partial 
differential equation in Wigner representation. For a derivation, we use the properties (|50"D 
and (|51"D of the operators T and H, to get from the master equation fl65|) for the density 
matrix elements p^ the corresponding operator equation 



+ 



N + l 



2fpf t - f f fp - pf f f ) + JV (2f t pf - f f f p - pff 1 ") } . 



(77) 



The dissipative part of this equation is the same as for the undriven dissipative harmonic 



oscillator |fL4|l , but with the shift operators for Floquet states instead of the usual creation 
and annihilation operators. 

Interestingly, the master equation in ([H]) now exhibits Lindblad form |T^JT5|1 . 

By substituting (|5"2"D, we get an operator equation which only consists of position and 
momentum operators. Transforming them into the Wigner representation, we find 



L(t) 



1 



pd x 



m 2 
with the coefficients 



(d x x + d p p) + k(t)xd x + 1 (D xx (t)d 2 x + D xp (t)d x d p + D pp (t)d 2 p ) (78) 



D xa! (t) = h£(t)£(t)(N + l/2)/m, 

D xp (t) = h (iimit) + &(t)Q(t)) (N + 1/2), 

D pp (t) = mh$>{t)i°(t)(N + l/2). 



(79) 
(80) 
(81) 



The fact that there are also dissipative terms in Eq. ([F8|) containing derivatives with 
respect to x is a consequence of the RWA: Its effect is equivalent to using instead of ([14]) the 
coupling Hamiltonian H^ A = J2v9v(abl + a^b^, where a and b v are the usual annihilation 
operators of the system and the bath, respectively. This introduces an additional coupling 
term oc pp u . In the next subsection we show how to avoid this RWA, by going back to the 
original Markov approximation, Eq. (p2|) . 



4- Fokker-Planck equation without rotating-wave approximation 

In the present case of a bilinear system, driven or not, for which the classical motion is 
integrable, the knowledge of the classical dynamics opens a more direct access also to the 
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quantal time evolution. Specifically, the interaction-picture position operator x(t,t') for the 
corresponding undamped quantum system is given by the solution of the classical equation 
of motion in the limit 7^0, indicated by the superscript °. In our case the classical solution 
is given by (|TT[) . The corresponding interaction-picture position operator reads 

i(t,f) = -x^^l + lG\t,t'). (82) 
at' m 

Inserting it into (0), we obtain a master equation in Markov approximation with respect to 
the quasienergy spectrum without expanding into Floquet states of the Schrodinger equation. 
Even with the rotating-wave approximation avoided, the resulting equation has already a 
simple structure: It is of the same form as the master equation derived in Sect. [II 1 B 1| , but 
with time-dependent transport coefficients 

roo roo 

7(t) = 2 7 / duu / dr sin(wr)G°(t - r,t), (83) 
Jo Jo 

(84) 



_ , . mh f°° , , ( %uj \ f°° . , , dG°(t-T,t' 

L>pp{t) = — — / dwcucoth [ — — J / dr cos(c<jr^ 



poo 

D 



'xp\ 




7r Jo V 2/crT J Jo dt' 

h 



t'=t 



ri f°° I huj \ r°° 

(*) = -/ dwwcoth — — / dr cos(wr)G°(t - r,t). (85) 
7r Jo \2k B T J Jo 

To evaluate these expressions, we substitute the undamped limit of Eq. (|T0|), 

G°(t, t') = c°c°, sin [fi°(t - t') + Q{nt - n't')] , (86) 

n,n' 

and exploit the sum rule (0) for the c n , to find, as in Sect. |III B 1| , 

j{t) = 7. (87) 

The explicit time dependence in G(t, t') results in a time dependence of the coefficients D pp 
and D xp . Averaging the transport coefficients over a period of driving, we find for D xp with 
the sum rule (0) again the expression (pT[), as in Sect. [Ill B 1| . Here, we have to choose the 
cutoff u>b much larger than the relevant frequencies /x° + nO. 
For D pp we find in an average over a period of driving 

Dpp= 1 -mh £ [cW + nn)]\oth H ^J n) . (88) 



Unlike the corresponding expression in the Sect. |111 B 1| , Eq. (p9|), the diffusion D pp now 
accounts explicitly for the quasienergies fr(/z° + nfl) instead of the energy huo. Thus the 
quasispectrum approach is reflected solely by a driving-induced modification of the momen- 
tum diffusion D pp . 

The Fokker-Planck equation for W(x,p, t) is now of the same structure as in the case of 
Markov approximation with respect to the unperturbed spectrum. Therefore the solution 
and the Floquet- Wigner functions remain the same, up to a different momentum diffusion 
Dpp. 

In contrast to the Fokker-Planck equation with RWA in the last subsection, the terms 
with d x x and d 2 x are now absent. In addition, the cross diffusion D xp in (R3) is completely 
different, and unrelated to the one in the RWA case fl8"0|). It originates from a principal part 
that has been neglected in the derivation of 
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D. Asymptotics 



1. The conservative limit 
In contrast to the Markov approximation with RWA in Sect. [Ill C 3|, the variances in 



both Markov approximations without RWA still depend on the friction 7. To obtain the 
conservative limit 7 — > of these, we insert the Green function (|TI]) into fl3"T|) and get 

7 £) / e -yt-i(2^+(n+n')Q)t 

0- xx (t) = -IT-f J2 °n C n' fl(t)' 



2m2 V 7 - i(2/x + in + n')n) 

■yt-i(n-n')flt p -yt+i(2fi+(n+n')n)t \ 

-2/i(t)/ 2 (t) 7 ^ +/ 2 2 (t) -7^ -7 - ,, n , . (89) 

7 — i(n — n'jil 7 + i(2/i + (n + n')\l) J 

In the limit of low damping, 7 |/i + nO| for any integer n, only the case n — n' of 
the second term in the brackets remains. Note that this condition is violated in parameter 
regions where the Floquet index becomes an multiple of Q, as is the case along the borderlines 
of the regions of stability in parameter space (cf. Fig. [1]). 
For the position variance, we get 

<r~(t)=A^&{t)${t), (90) 

where 

00 9 

A= E (c°n) (91) 

n=— 00 

denotes a number of order unity. 
In an analogous way, we find 

a xp (t) = + Cx ^)) , (92) 



a 



pp 



(t) = AD pp £° 1 (t)g(t). (93) 



Besides the prefactor, these variances are the same as for the master equation with RWA in 
Sect. II7TC1 

Moreover, in this limit 7^0, all diagonal elements W nn (x,p,t) are Floquet functions 
with the quasienergies \x nn = 0. However, they are different from the Wigner representation 
of the stationary solutions flT3"| ) of the corresponding Schrodinger equation, which are of 
course solutions of ( |33"D with 7 = 0. Due to the degeneracy of the Floquet indices, this is no 
contradiction. The W nn (x,p,t) can be viewed as dissipation-adapted Floquet functions. 

For consistency, we check the uncertainty relations for the asymptotic solution. It is 
satisfied if the variances fulfill the inequality 



0~xx 

(t) a xp (t) 

Vxp{t) <Tpp(t) 



fD pp A\ 2 2 



V m 

which we have verified numerically for the case of the Mathieu oscillator 



> (94) 
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2. The high-temperature limit 



In the limit of high temperatures k B T S> hoo-o, we expect the Fokker-Planck equation for 



the Wigner function to give the Kramers equation for the classical Brownian motion |p28 
i.e. an equation of the form ([33|) with the diffusion constants D xp = and D pp = mk^I '. 

In the standard approach (Sect, [ill B| ) and the quasispectrum approach without RWA 
(Sect. |1HC4|) , the Fokker-Planck equation is already of the required structure. With = 
C pT| the cross diffusion D xp vanishes in the high-temperature limit. For D pp , we use 
cothx = 1/x + 0(x) and get 

D pp = mk B Tj2(c°n) 2 (v° + nfi). (95) 

n 

With the sum rule (0), this reduces to D pp = mk B T. 

In the quasispectrum approach with RWA in Sect. [Ill C 3| , the variances and diffusion 



constants scale with N + 1/2. This factor reads, in the high-temperature limit, 

^ + i = ?((S )^ = ^. ,96) 

Therefore the diffusion constants D xx and D xp remain finite and the Fokker-Planck operator 
( f78l) does not approach the Kramers limit for high temperatures. Nevertheless the asymp- 
totic variances in RWA coincide for high temperatures, with the classical result in the limit 
7 -> 0. 



IV. NUMERICAL RESULTS 

In this section, we compare our approximate results to exact ones, obtained from the 
path- integral solution in Ref . [|IJ . Specifically, we give the numerical results for the Mathieu 
oscillator, i.e., we use 

k{t) = m (wq + e cos Vtt) . (97) 
This is an experimentally important case in view of the fact that it describes the Paul trap 

By inserting ( |9^) and the ansatz (|5]) into (^), we obtain the tridiagonal recurrence relation 
£ Cn _! + 2 (uol - 7 2 /4 - (/i + nnf) c n + ec n+1 = 0. (98) 



From this equation, the classical Floquet index fi and the Fourier coefficients c n are deter- 
mined numerically by continued fractions |[24|| . 

In the figures we use the scaled quantities i = Qt/2, u Q = 2u /Q and e = 2e/Q 2 . The 
external period thus takes the value T = ir. Position and momentum are scaled via x = 
(2h/mVL)~ l l 2 x and p = (mhQ/2)~ 1 / 2 p, respectively. The overbar for the scaled quantities 
has been suppressed in the figures. 

The influence of the quasienergies on the equation of motion fl33|) is given by differ- 
ent diffusion coefficients D pp . In Fig. |2], we compare the momentum-diffusion coefficients 
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between the Markov approximation with respect to the unperturbed spectrum, given by 
Eq. (|29|), and the Markov approximation that relates to the quasienergy spectrum, given by 
Eq. (55). We have scaled the values to the classical momentum-diffusion coefficient ink^T . 
The parameters u 2 . and e are varied along the full line in the inset. Note that within the 
unstable regimes, perturbation theory is not valid. Nevertheless, Eq. flSSJ) gives a smooth 



interpolation. The discrepances become most significant for strong driving and large uJq. 
For both, low driving amplitude eCuJ and high temperature T ^> fruo/kB, the difference 
vanishes. 

The variances u xx {t) and cr pp (t) of the Markov approximations without RWA are com- 
pared against the exact results jl| in the panels |3]a and [5]b. The chosen driving parameters 



2 = 6.5 Q 2 and e = 7Q 2 lie inside the fifth stable zone (/i = 4.53513 Q/2). The temperature 
k^T = 0.5 Ml is sufficiently large, but with quantum effects still appreciable. We note that 
the improved Markov treatment in Sect. |III C 4| , that accounts for the quasienergy differ- 
ences, agrees better with the exact prediction. In the Figure we depict asymptotic times 
t > 50, where transient effects have already decayed. The asymptotic co- variance elements 
retain the periodicity T = ir of the external driving. For the chosen parameters the relative 
error is reduced by the use of the improved Markov scheme by approximately 30%. 

The relative error r) xx (t) = cr^ arkov (i) / a™ ct (t) of the position variance for these two 
Markov approximations is depicted in panel [|. Note that the maximal deviations do not 
occur in the extrema, but happen to occur in the regions with negative slope. 

As depicted with Fig. [5], the quality of both Markov approximations worsens with in- 
creasing dissipation strength 7. This reflects the breakdown of the weak coupling approach 
when strong friction is ruling the system dynamics. 

Results for the Markovian treatment within RWA, given in Sect. Ill C 3 , are depicted for 
the position variance <J xx (t) in Fig. The driving parameters are the same as in Fig. |3|. 
For this example, the quality of agreement to the exact result is similar for both Markov 
approximations. Nevertheless, the solution without RWA yields — up to a scale — a better 
overall agreement with the exact behaviour over a full driving period T. 



V. CONCLUSION 



We have used the parametrically driven harmonic oscillator as a simple working example 
to compare various versions of the Markovian approach to the quantum dynamics of period- 
ically driven systems with dissipation, and to provide a synopsis of a number of alternative 
representations, each of which emphasizes different aspects of the same underlying physics. 

The principal distinction to be made among possible Markovian approaches to the driven 
dissipative dynamics, refers to the degree to which changes in dynamical and spectral proper- 
ties of the central system due to the driving are taken into account. In the crudest treatment, 
the non-unitary terms in the master equation are derived ignoring the explicit time depen- 
dence of the Hamiltonian, and the driving appears only in the unitary term. An improved 
master equation is obtained if the central system and the driving are coupled to the heat 
bath as one whole. The energy-domain quantity relevant for all subsequent developments 
is then the quasienergy spectrum, obtained within the Floquet formalism, instead of the 
unperturbed spectrum. In the time domain, the quantities entering the dissipative terms 
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of the master equation, such as Heisenberg-picture operators of the central system, gain an 
explicit time dependence with the periodicity of the driving. As a bonus, the Floquet treat- 
ment of the central system with driving yields a well-adapted basis, the set of eigenstates 
of the Floquet operator. Representing the master equation in this basis completely removes 
the unitary term. 

Besides the differences in representation, the use of the improved Floquet-Markov ap- 
proximation in Sect. [Ill C 4| results mainly in a modified momentum diffusion that dependss 
on the quasienergy spectrum instead of the unperturbed spectrum of the central system. 
The difference becomes significant in the limits of strong driving amplitude and low temper- 
ature. An additive time-dependent external force, applied in addition to or instead of the 
parametric driving, undergoes a renormalization which vanishes, however, in the case of an 
Ohmic bath. 

Even within the improved Markov approach, finer levels of approximation can be distin- 
guished. A significant simplification of the master equation is achieved by a rotating-wave 
approximation, i.e. here, by neglecting reservoir-induced virtual transitions between Floquet 
states of the central system that violate quasienergy conservation. The resulting master 
equation has Lindblad form, with creation and annihilation operators acting on Floquet 
states, and thus manifestly generates a dynamical semigroup. This is not the case if the 
RWA is avoided. Apparently a drawback, the lack of a Lindblad structure in the master 
equation without RWA faithfully reflects the failure of the Markov approximation on short 
time scales. 

An analogous situation as with the Lindblad form of the master equation arises with 
its Floquet structure. If all coefficients are at most periodically time dependent, then the 
equation of motion for the reduced density operator complies with the conditions for ap- 
plicability of the Floquet theorem. As a consequence, the solutions can be cast in Floquet 
form, i.e., can be written as eigenfunctions of a generalized non- unitary Floquet operator 
that generates the evolution of the density operator over a single period. Since all variants 
of the Markov approximation discussed herein truncate the memory of the central system on 
time scales shorter than the period of the driving, the corresponding master equations have 
Floquet structure throughout. The exact path- integral solution, in contrast, allows for mem- 
ory effects of unlimited duration and thereby generally prevents the consistent definition of 
a propagator over a single period. 

Additional insight is gained by discussing the dynamics in terms of phase-space distri- 
butions, specifically, in terms of the Wigner representation of the density operator and its 
equation of motion. In this representation, the Floquet formalism is a useful device to con- 
struct and classify solutions. Since all Fokker-Planck equations obtained are time periodic, 
as are the corresponding master equations, their solutions may be written as eigenstates 
of a Wigner- Floquet operator (the Fokker-Planck operator evolving the Wigner function, 
integrated over a single period), or Wigner- Floquet states in short. They represent the 
quasiprobability distributions closest to the Floquet solutions of the corresponding classical 
Fokker-Planck equation. 

Wigner-Floquet states with a purely real quasienergy correspond to asymptotic solutions. 
They are not literally stationary but retain the periodic time dependence of the driving. 
Since we are here dealing with a linear system, the asymptotic quasiprobability distributions 
follow the corresponding classical limit cycles. In the case of parametric driving, these limit 
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cycles are trivial and correspond to a fixed point at the origin. A time dependence arises 
only by the periodic variation of the shape of the asymptotic distributions. 

Concluding from a numerical comparison of certain dynamical quantities, for the specific 
case of the Mathieu oscillator, the attributes "simple" and "improved" for the two basic 
Markovian approaches prove adequate. Results for the Markov approximation based on 
the quasienergy spectrum show consistently better agreement with the exact path-integral 
solution than those for the Markov approximation with respect to the unperturbed spectrum. 
However, even in parameter regimes where the respective approximations are expected to 
become problematic, the differences in quality are not huge and the agreement with the exact 
solution is generally good. Technical advantages of the Markov approximation in general 
and of its various ramifications — easy analytical and numerical tractability, desirable formal 
properties such as Floquet or Lindblad form of the master equation — can justify to accept 
their quantitative inaccuracy. 
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APPENDIX A: SOLUTION OF THE CHARACTERISTIC EQUATIONS 

In this appendix, we solve the equation of motion for the Wigner function by the method 
of characteristics. For simplicity, we use here units with m — 1. We write W(x,p, t) as 

W(x,p,t) = J dXdPe^+^e 5 ^'^. (Al) 

By this ansatz, equation ( |3"B| ) is transformed to the quasilinear partial differential equation 

F(X,Sx,P,S P ,t,S t ) = (A2) 

for S(X, P, t), where T is given by 

T=S t - XS P + jPSp + k{t)PS x + jD pp P 2 + ~?D xp XP. (A3) 

We denote the partial derivatives of S(X,P,i) with respect to X, P and t by Sx, Sp and 
S t , respectively. 



The characteristic equations |26] of (|A2|) are given by 



*=f£ = l> ( A4 ) 



dS t 
?_ 

dS x 



X = ^ = k(t)P, (A5) 
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p= w r =lP - x ' (A6 » 



dX 



Sx = -— = S P - 1 D xp P, (A7) 



S P = = -iSp - k(t)S x - ^D^P - 1 D xp X } (A8) 

* = -£ = -^ < A9 > 

whose solutions give the characteristics of the partial differential equation ( |A2| ). 

Equation ( |A4] ) means that the characteristics can be parameterized by the time t. Instead 
of equation (|A9|) , we will use (|A^) to get an expression for S t . So we only have to solve 
(|A5|)- (|A8|) . The solutions of these equations can be traced back to the fundamental solutions 
fi(t) of the classical equation of motion @). 

From ( fA5| ) and ( |A6| ) , we find 

P - 7P + k(t)P = 0. (A10) 

This is simply the classical equation of motion with a negative damping constant. Therefore 
the solutions for X and P read 

P(t) = -c 1+ e*f 2 (t) + C2+e^/i(t), (All) 
X{t) = c 1+ e*/ 2 (f) - ca+e*/^), (A12) 

where Cj + denote integration constants. 
From (|A7|) and (TAB]) we find for 5* x 

£x + jS x + k(t)S x = -2 7 DP, (A13) 

which is the classical equation of motion with an inhomogeneity. The effective diffusion 
constant D is given by 

D = D PP + 1 D xp . (A14) 
With the integration constants q_, we integrate ( |A13| ) with the Green function (||) to 

Sxit) = ci_/i(t) + P2-/2W - 2 7 P f &t'G(t,t')P(t'), (A15) 

J to 

and get by use of ( |A7|) 

5 P (t) = Cl ^j\(t) + c 2 J 2 (t)-2 1 Dj t t dt'^^P(t')+ 1 D xp P(t). (A16) 



By inserting 



P(t') = G(t, t')X(t) + ^f^P(t), (A17) 
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obtained from Eqs. ( |A11|) and ( |A12|) , we get a result for Sx and Sp that depends only on 
the endpoints of the characteristics. Now together with Eq. (A2), we have an expression for 
gradS^X, P, t) = (Sx, Sp, S t ), which can be integrated to 



with 



S(X,P,t) = (ci_/i(t) + c*-h(t))x + (ci_A(t) + C2_/ 2 (t))P 
— n a xx(t, to)X 2 — cr xp (t, to)XP — -<J pp (t, to)P 2 , 



a xx (t,t ) = 2 1 D f dt' [G(t,t')f, 

J to 

a xp (t,t ) = 2 1 D J* dt'G(t,t')^-G(t,t') 



Vp P {t, t ) = -~fD xp + 27 D [ dt' 

J t 



d_ 

dt 



G(t,t') 



(A18) 

(A19) 
(A20) 

(A21) 



By inserting S(X,P,t) into (|A1|) , we find a solution for the Wigner function W(x,p,t). 

The integration constants q± are of course constant along the characteristics. Therefore 
the Poisson brackets between the expressions Ci±(X, Sx, P, Sp, t) and T[X, Sx, P, Sp, t, S t ) 
vanish By transforming back from Fourier space to real space, one finds that the 

operators q± = Ci±(—id x ,—ix,—id p ,—ip,t) commute with the operator d t — L(t), whose 
nullspace is the solution of the equation of motion. Therefore, the q± are shift operators in 
the subspace of solutions, i.e. if W(x,p,t) is a solution of fl33|), then &i±W(x,p,t) is also a 
solution. 

For the q± we find 



C 2 - 



\{hit)d x + hm 

Uf2(t)d x + 



(t,t )d x + a xp (t,t )d p 
-if 2 (t) (p + a xp (t, t )d x + a pp (t, t )d p 
-ifi(t)(x + a xx (t,t )d x + a xp (t,t )d p 
+ifi(t)(p + a xp (t,to)d x + a pp (t,to)dj 



(A22) 
(A23) 

(A24) 
(A25) 



Note that because of the linear structure of the characteristic equations, there is no 
ambiguity concerning the ordering of operators. 

The operators Q i+ (t), used above, are proportional to the c i+ . 



APPENDIX B: THE ADDITIVELY DRIVEN HARMONIC OSCILLATOR 

In this appendix we present the Markovian master equation within the quasi-spectrum 
approach when the parametric oscillator is subjected to additional additive driving —xF(t), 
i.e. 
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H(t) = H s (t)-xF{t). (Bl) 

With Hs{t) being a time-independent harmonic oscillator, i.e., k(t) = tjkJq, the correspond- 
ing Markovian master equation in RWA for the dissipative system has already been given in 
H. Herein we generalize these results for the combined time-dependent system Hamiltonian 



in dBT|) 



It is known that the only effect of the driving force F(t) on the (quasi-) energy spectrum 
of a parametrically driven harmonic oscillator is an overall level shift |§. Thus the level 
separations remain unaffected and we expect no change in the dissipative part of the master 



equation (|27|) . 

The classical equation of motion, which is also obeyed by the interaction-picture position 
operator, now reads 

mx + k(t)x = F(t), (B2) 
and can be integrated to yield the interaction-picture position operator 

x(t,lf) = _a gG y) + P-G°(t,t') + - f dt" G°(t,t")F(t"). (B3) 
ot m m Jt' 

Thus we obtain a c- number correction to the interaction-picture position operator (|32"1) , given 
by the third term. After inserting (p3|) into (0), the generalized Markov approximation 
emerges as 

PS = ... + j;F(t)[x,P8] (B4) 

-ilE^/ dTA„(T)[x, Ps \- / dt'G (t-r,t')F(t'). (B5) 
n v Jo m Jt 

The dots denote the old result for F(t) = 0, given by the right hand side of Eq. (^7]). The 
term in the first line stems from the reversible part of the master equation (p2|); the second 
one is a correction of the driving force due to the interaction with the bath. Thus the 
equation of motion for the density operator has the structure 

P8 = ... + ^F(t)[x,p8] (B6) 
with an effective total driving force 

F(t) = F(t) + / dul(u) / dr sincjr / dt'G°(t - T,t')F(t'). (B7) 

rrm Jo Jo Jt 



Note that the dissipative parts of ( |B6|) are not affected by the additive driving force F(t). 
This makes explicit, that we must use a parametric time-dependence to study differences in 
the dissipative parts resulting from the Markov approximation with respect to the energy 
spectrum versus the Markov approximation with respect to the quasienergy spectrum. 

With an Ohmic bath, I{uj) = rwyu, the integral in ( |B7|) vanishes and we obtain F(t) = 
F(t). Thus in contrast to an explicit parametric time dependence k(t) in the quadratic part 
of the Hamiltonian, the time dependence of an additive force, in this case, does not change 
the Markovian master equation of the dissipative system. 
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FIGURES 




FIG. 1. Stability of equation @ with 7 = for the case of a Mathieu oscillator. In the white 
areas the Floquet index fi is real, which corresponds to stable solutions. In the shaded areas fi 
is complex and therefore one of the fundamental solutions (||) is unstable. On the borderlines \x 
becomes a multiple of Q./2. 




FIG. 2. The diffusion constant D pp for the simple (dotted) and the improved (dashed) Markov 
approximation compared to the time-average of the exact value in units of the classical diffusion 
constant mk-QT for k^T = 0.5 HQ. The parameters cjq and e are indicated by the full line in the 
insert (units as in Fig. [l]). 
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FIG. 3. The asymptotic variances cr xx (t) (a) and cr pp (t) (b) with period T = lixjQ, for the 
simple (dotted) and the improved (dashed) Markov approximation, compared to the exact result 
(full line) for the parameters e = 7Q 2 , ujq = 6.5 f2 2 , k^T = 0.5 and 7 = 17/20. The scaled 
driving period T = ir is indicated in panel (a). 
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S 1.3 




FIG. 4. Relative error rj xx (t) = o'^, arkov (t)/cr^ act (i) for the position variances of panel |a. 




FIG. 5. The time averaged variance cr xx (t) for the simple (dotted) and the improved (dashed) 
Markov approximation, compared to the exact result (full line) for the parameters e = 
= 6.5 n 2 and k B T = 0.5 ftO. 
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FIG. 6. Position variances obtained with the Markov approximation with respect to the 
quasienergy spectrum with (dotted) and without (dashed) RWA, compared to the exact result 
(full line) for 7 = 17/10 for k-gT = 0.5 Ml The driving parameters are e = 7Q 2 and oJq = 6.5 . 
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